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We study the superconducting phase with two component order parameter scenario, such as, 
dx2_y2 + e' Sa, where a = xy, x + y . We show, that in absence of orthorhombocity, the usual 
dx2_y2 does not mix with usual Sx2+y2 symmetry gap in an anisotropic band structure. But the Sxy 
symmetry does mix with the usual d-wave for 9 = 0. The d-wave symmetry with higher harmonics 
present in it also mixes with higher order extended s wave symmetry. The required pair potential to 
obtain higher anisotropic d^2_y2 and extended s-wave symmetries, is derived by considering longer 
ranged two-body attractive potential in the spirit of tight binding lattice. We demonstrate that 
the dominant pairing symmetry changes drastically from d to s like as the attractive pair potential 
is obtained from longer ranged interaction. More specifically, a typical length scale of interaction 
^, which could be even/odd multiples of lattice spacing leads to predominant s/d wave symmetry. 
The role of long range interaction on pairing symmetry has further been emphasized by studying 
the typical interplay in the temperature dependencies of these higher order d and s wave pairing 
symmetries. 
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I. INTRODUCTION 



Many experiments were performed to find clues regard- 
ing mechanism of high- Tc superconductivity and the na- 
ture of the superconducting pair wave function. Notwith- 
standing this effort the nature of the orbital symmetry 
of the order parameter is not yet known completely af- 
ter a decayed of its discovery although strong evidence of 
a major dx2_y2 symmetry exists [ll-p|. Phase and node 
sensitive experiments also reported a sign reversal of the 
order parameter supporting d wave symmetry [H. The 
most current scenario as appears fronr various experi- 
ments and theory that the pairing symmetry of these 
family could be a mixed one like dx2_y2 -|- e'^a where a 
could be something in the s wave family or d^y The elec- 
tron doped Nd2-xCexCu04 superconductors are however 
pure s wave like H . 

Tunneling experiments had questioned the pure d-wave 
symmetry M as the data were interpreted as an admix- 
ture of d and s-wave components due to orthorhombic- 
ity in YBCO WM- Possibility of a minor but finite 



idxy symmetry alongwith the predominant dr^2_y2 has 
also been suggested pO[ | in connection with magnetic de- 
fects or small fractions of a flux quantum $o = hc/2e 
in YBCO powders. Similar proposals came from various 
other authors in the context of magnetic field, magnetic 
impurity, interface effect etc. |ll^-0] These proposals got 
the correct momentum when experimental data on lon- 
gitudinal thermal conductivity by Krishana et al, |l2] 
of Bi2Sr2CaCu208 compounds and that by Movshovich 
et al, |l3] showed supportive indication to such propos- 
als. There are experimental results related to interface 
effects as well as in the bulk that indicates mixed pairing 



symmetry (with dominant d-wave) p3{ , thus providing a 
strong threat to the pure d wave models. 

In this paper our main aim is to study the possibil- 
ity of a mixed pairing symmetry state with A(fc) — 
Ad 2_ 2 + e'^Sa where a = xy,x^ + y^ for 9 = 0, 7r/2 
with both d and s on an equal footing. We show that 
dj.2_y2 can mix with s^y in the tetragonal group for 9 — Q 
but not for 9 — 7r/2. The phase of the second conden- 
sate state is thus extremely important. We then show 
that eventhough the lowest order dj.2 _y2 cannot mix with 
.Sx2^y2, the corresponding higher order symmetries can 
mix freely with each other. By lowest order we mean the 
usual d-wave {i.e, simple coskx — cos ky form), extended 
s wave (i.e, simple cosfc^, -f cos ky form) and so on. By 
higher order we mean such symmetries with higher har- 
monics present in it, like cos^k^ ± cos ^ky form where 
^ = na (n — 1,2,3..) or even more complicated like 
cos 2kx cos ky ± cos kx cos 2ky and so on. This will be 
clearer as we proceed. Now, in order to obtain such pair- 
ing symmetry in the respective channels one needs effec- 
tive attractive pairing potential V{^k,£^k'). We derive, 
in the spirit of tight binding longer range attraction than 
the usual nearest (or next nearest) neighbour one such 
interaction potential. The potential V{(,k,£_k') therefore, 
changes the position of its minimum from that of the 
usual d or s wave cases for n > 1. We show, depending 
on the position of the pair potential or in other words, 
longer ranged attractions ^ = 2a, 3a, 4a etc. the domi- 
nant symmetry changes from dx^._y2 for ^ = a to s like 
otherwise. 

This study can particularly be justified based on the 
following grounds, (i) On general grounds, long range 
interaction arise from a decrease in screening as one ap- 
proaches the insulator. In specific models of supercon- 



ductivity like the spin-fluctuation mediated models, an 
increase in the antiferromagnetic correlation length oc- 
curs with underdoping. (ii) One of the potential theo- 
ries of high temperature superconductivity that favors d 
wave symmetry is the spin fluctuation theory p3. The 
gap symmetry of the spin fluctuation theory is however 
not the simplest d-wave but higher order d-wave, approx- 
imately of the form (cosfc^^ — cosfcj^)(cosfca; + cos ky)^ 
[ p5[ . Explicit fc-anisotrpy of the gap in spin fluctuation 
mediated superconductivity was obtained by Lenck and 
Carbotte [|l^,pl in BCS theory with the phenomenolog- 
ical spin susceptibility as pairing interaction using fast- 
Fourier- Transform technique, without any prior assump- 
tion about the symmetry of the gap. They concluded, 
the gap although have nodal lines along k^ — ky, does 
not have the simplest d-wave symmetry but rather higher 
order d wave symmetry with higher harmonics present in 
it. Therfore, this work provide a real space derivation of 
a pair potential that produces higher order d-wave sym- 
metry similar to that present in the spin fluctuation the- 
ory. (Hi) In the magnetic scenario of the cuprates [[l7| , 
one can set ^ equal to the magnetic coherence length 
which is larger than the lattice spacing |3|. The coher- 
ence length in the superconducting state which is differ- 
ent for different materials may be because a short range 
interaction requires larger densities than a long range one 
in order to produce coherent motion that leads to super- 
conductivity. (The Tc — X relationship is not unique in 
all high Tc systems, some starts to superconduct with 
very small doping, x whereas some systems require larger 
x). (iv) The high Tc systems are in very complicated 
circuit and the electronic correlation effects may not be 
adequately accounted unless one considers next nearest 
or further neighbour repulsion. Therefore, in the spirit 
of tight binding lattice the effective attraction may only 
arise with more distant attractive interation. (w) In a 
most recent angel resolved photoemission (ARPES) ex- 
periment by a well known group |18| , such requirement of 
long range interaction was realized. One of their essential 
findings is, as the doping decreases, the maximum gap in- 
creases, but the slope of the gap near the nodes decreases. 
This particular feature although consistent with d wave 
but cannot be fit by simple cos(2(/)) but requires a finite 
mixing of cosidcj)) as well, where is angle between k^, ky 
given as, tan'^ (ky/kx). The cos(60) contains higher har- 
monics than simple (cos k^ — cos ky). Rest of the lay out 
of the paper is as follows. In section II, we derive the pair 
potential required for higher anisotropic d and extended s 
wave symmetries. We also provide a brief prescription of 
finding coupled gap equations for the amplitudes of such 
higher anisotropic symmetries. In section III, we present 
and discuss in details all the numerical results providing 
strong signature of change in dominant pairing symmetry 
with range of interation. Finally, we conlude in section 
IV. 



II. MODEL CALCULATION 

Let us consider that the overlap of orbitals in differ- 
ent unit cells is small compared to the diagonal overlap. 
Then in the spirit of tight binding lattice description, the 
matrix element of the pair potential may be obtained as, 



V{q) 



^ F^e*?«^ = Vo^ + V,f\k)f{k') + V^g{k)g{k') 

5 

+F2/'*- {k)f^-y [k') + v^2r- (fc)r- (fc') 

+V:if{2k)f{2k') + F35(2fc)g(2fc') 
+2VJtm~ft{2k') + 2VJi{2k)~fi{2k') 
+2V4gi{2k)gi{2k') + 21/4.92(2^)32 (2fc') 
+V5f'^''^i2k)f^^{2k') + Vr,f''^{2k)f''-y{2k') 
+Vef''{3k)f{3k') + Veg{3k)g{3k') (1) 



where in the first result of the equation (nh Rs locates 
nearest neighbour and further neighbours, 5 labels and 
V„, n = 1, ...,6 represents strength of attraction between 
the respective neighbour interaction. The first term in 
the above equation Vq refers to the on-site interaction 
which is considered as repulsive but can be attractive as 
well giving rise isotropic s wave. In this paper, we shall 
not consider the isotropic s wave for a mixed symmetry 
with d wave (cf. py]). The form factors of the potential 
are obtained as, 

/ (nk) — cos{nkxa) — cos{nkya) 

g{nk) — cos{nkj:a) + cos{nkya) 

f ^^{nk) — 2sm(nkxa) sm(nkya) 

f'^^^{nk) — 2 cos{nkxa) cos{nkya) 

fi{2k) = cos{2kxa) cos{kya) — cos(fc2;a) cos{2kya) (2) 

/2 (2fc) = sm(2kxa) sin(fcya) — sin(fca;a) sin(2fcj,a) 

gi{2k) ~ cos{2kxa) cos(fcj,a) + cos(fca;a) cos(2fcj,a) 

g2{2k) = sm{2kxa) sm{kya) + sin(fca;a) sin(2fcj,a) 

where f'^{nk), g{nk) leads to usual dx^-y^, Sx^+y'^ pairing 
symmetry for n — 1 and unusual or higher order dj,2_j^2, 
Sx^+y^ pairing symmetry respectively which results from 
interations along the x and y axes {i.e, f**, 3'"'^, 6*'* neigh- 
bour interaction). While the usual and higher order dxy, 
Sxy pairing symmetry results from f'^''y{nk), f^'^^ink), 
the 4*'' neighbour interaction gives rise to unconventional 
d and extended s-wave pairing symmetry through f^(2k) 
and gn{2k) given in equation (2). In deriving Eqs. (|l|,2) 
terms responsible for triplet pairing which are not im- 
portant for high Tc systems are neglected. We shall 
discuss now the mixed phase symmetry of dx2^y2 with 
other symmetries taking two of the potential terms at 
a time, namely, a combination of potential terms in (Q) 
(2"''^, S'''*), (6*'', 7*'*), (14*'', 15*'') gives rise to pairing sym- 
metry A(fc) = Ad^,_^Mf'm + e''^s^.^,A0)9m 



where ^ = na, a is the lattice constant and will be 
taken as unity. Similarly, a comibnation of (2"'^,4*''), 
(6*^, 12*'') and so on will give rise to pairing symmetry 
A(fc) = Arf^,^^,(0)r(efc) +e*^Arf^„(0)r-(efc) etc. 

Free energy of a superconductor with arbritary pairing 
symmetry may be written as, 



III. RESULTS AND DISCUSSIONS 



F,,,, = -^ E i^^(i 



/3 



+ e 



-pPEk 



A, 



/c,p— ± 



Vkk' 



(3) 



where Ek — y/i^k — m)^+ I ^fe P are the energy eigen 
values of a Hamiltonian that describes superconductiv- 
ity. We minimize the free energy, Eq. (H) i.e, dF/d \ A | 
= 0, to get the gap equation as, 



Afe = E Vkk' 



2Ek 



■tanh(^— ) 



(4) 



where e^ is the dispersion relation taken from the ARPES 
data ||9(] and /i the chemical potential will control band 
filling through a number conserving equation given be- 
low. For two component order parameter symmetries 
as mentioned above, we substitute the required form of 
the potential and the corresponding gap structure into 
the either side of Eq. (|j) which gives us an identity 
equation. Then separating the real and imaginary parts 
together with comparing the momentum dependences on 
either side of it we get gap equations for the amplitudes 
in different channels as, 



T^y. 



ML 

2Ek 



tanh 



pEk 



1,2 



(5) 



Considering mixed symmetry of the form A(/s) = 
^d,2_„2 iO)f'^{nk) + As^2^^2 {0)g{nk) one identifies Ai = 
a/,_[,(0), A2 = aJ^7^(0) and fl = f'ink), fl = 
g{nk). Similarly, for mixed symmetries of the form 
A(fc) = Ad , , {0)f{nk) + A,^„ (0)/"- where a = s,d 

A2 — Aq^j^(O) and /| = j"^" and so on. The potential 
required to get such pairing symmetries are discussed in 
Eq. ®. 

The number conserving equation that controls the 
band filling through chemical potential, fi is given by. 



k 



1- 



{ek 



Ek 



^)tanh^ 



(6) 



We solve self-consistently the above three equations 
(Eq.H and Eq.@) in order to study the phase diagram of a 
mixed order parameter superconducting phase. The nu- 
merical results obtained for the gap amplitudes through 
Eqs. (0Jq) will be compared with free energy minimiza- 
tions via Eq. (H) to get the phase diagrams. 



We present in this section our numerical results for a 
set of fixed parameters, e.g, a cut-off energy nc= 500 K 
around the Fermi level above which superconducting con- 
densate does not exist, a fixed ratio V1/V2 = 0.71 in Eq. 
(o) between the strengths of pairing interaction channels 
through out. In figures 1 and 2 we present results for 



Aik)^A,,JO)f'^m+e^<^A, 



, {0)g{^k) symme- 



tries for = 7r/2 and 9 — respectively. Such symmetries 
would arise from a combination of two component pair 
potentials (2"'', 3''''), (6*'', 7*''), (14*'*, 15*'') and so on. 
We shall discuss only the results oi 9 = and 9 = tt/2. 
These two phases of 9 can cause important differences (cf. 
figures 3, 4). It is known that for any 9^0, time rever- 
sal symmetry is locally broken |19[| which correspond to a 
phase transition to an almost fully gapped phase (except 
at the points ±7r/2, ±7r/2 due to common nodal points 
from both the channels) from a partially ungapped phase 



of d,r 



symmetry. On the other hand, the ^ = phase 



still remains nodeful, although the nodal lines shifts a lot 



from the usual k^. 



ky lines of the d^2_ 
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FIG. 1. Amplitudes of the Ad , , (solid lines) and 
As 2 2 (dashed lines) as a function of band filling p for 
6 = 7r/2 (i.e, A^2_y-2 + is3j2_|_y2) phase in various values of 
f/a. While the usual A^2_y-2 does not mix with usual Sx2+y2 
(a), higher component dx2_y2 and 8^2 +y2 (c), (d) can mix 
with each other freely even in absence of orthorhombocity. 
It is worth noticing that the change in the dominant pair- 
ing symmetry with ^/a {e.g, for ^/a = 2 the only dominant 
symmetry is s wave like). 



FIG. 2. Same as that of figure 1 except 6 = Q {i.e, 
dx2_y2 -|-Sx2^y2 symmetry). The predominant symmetry al- 
ways tries to expel (minimize) occurance of the other symme- 
try at its optimum doping. 



The solid lines represent the amplitude oid^2_y2 chan- 
nel whereas the dashed lines indicate that of Sx2^y2. 
These figures (1 & 2) clearly demostrate that the usual 
dx2-y2 and Sr^2_^_y2 symmetries do not mix with each other 
(cf. Figures 1(a), 2(a)) but the higher order d^2_y2, 
s^2_^_y2 symmetries do mix with each other (cf. figures 
l(c,d), 2(c,d)). In fact, as the interaction becomes longer 
ranged (i.e, ^/a — 1,2,3,4 as is demonstrated in figures 
1, 2 (a), (b), (c), (d) respectively) the dominant sym- 
metry changes drastically; as the typical length ^ is odd 
multiple of the lattice constant, the dominant symmetry 
at lower doping is dx2_y2 like whereas when the ^ is even 
multiple of the lattice constant, the dominant symmetry 
at lower doping is something in the s-wave family (see 
also figures 3, 4). 




0.010 h 

0.008 

0.006 

0.004 h 

0.002 

0.000 



As the typical length ^ is increased the predominant 
symmetry at the optimal doping CQ]) changes from d- 
wave at 5 = a, to an extended s-wave Sx2^y2, Sxy for 
^ = 2a, to again a predominant d-wave symmetry at 
^ = 3a and finally for ^ = 4a to extended s wave symme- 
try for 9 = 7r/2. These phase diagrams (figures 1,2,3,4) 
drawn at T = 1 mK does not change the scenario even 
for = 0, in the mixed phase of d-wave with s^2_|_j,2 
symmetry but causes significant change for that with Sxy 
symmetry (cf. Fig. 4). More significantly, the case of 
^ = 2a is universal (i.e, independent of 9 and s^2_|_j,2 
or Sxy mixing with d -wave), the dominant symmetry 
at zero temperature is s-wave type. This work there- 
fore, has revealed in a significant way the change in pre- 
dominant pairing symmetry as the interaction range is 
changed at T=0. It is to be noted that in contrast to 
hole doped material, the electron doped materials (like 
Nd2-xCexCu04) have no signature of dominant d-wave 
symmetry. Furthermore, the antiferromagnetic phase in 
the electron doped systems is more extended or exists till 
larger doping in comparison to the hole doped material. 
Therefore, considering models related to spin fluctuation 
mediated superconductivity, the longer range attraction 
should be more important. In the present picture, we 
showed that such longer range interaction cause change 
in the pairing symmetry which might make this study to 
have imporatant bearings for the high-Tc compunds. 
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FIG. 3. Amplitudes of the Ad 2_ 2 (solid lines) and As^y 
(dashed lines) as a function of band filling p for 8 — -k 12 
{i.e, dx2_y2 + iSx2+y2) phase in various values of ^/a. While 
the usual dx2_y2 does not mix with usual Sxy (a), higher 
anisotropic dx2_y2 and Sxy (c), (d) can mix with each other 
freely even in absence of orthorhombocity. It is worth notic- 
ing that the change in the dominant pairing symmetry with 
the typical length ^/a {e.g, for ^/a — 2 the only dominant 
symmetry is s wave like). The figure (a) should particularly 
be contrasted with that of figure 4. 



Some interesting features of the data presented is that 
optimal doping remains unchanged irrespective of ^ that 
causes a significant crossover in the dominant symmetry 
of the order parameter. The position of the d-wave does 
not change appreciably except the case of ^/a = 4 while 
the extended s wave region moves drastically with ^. In 
particular, for ^/a ~ 1, the extended s wave family has 
finite amplitude only at densities close to zero {p ~ 0) (cf. 
figures 1,2,3) leading to no mixed phases except the out- 
standing case oi 6 = for Sxy (cf. figure 4). In ^/a — 2 
case, the extended s-wave family completely takes over 
the position of the d-wave that it had in case of ^/a = 1. 
For ^/a — 3, the d-wave regains its poisition although 
both the amplitude and width decreases to about 50% to 
that of the S,/a = 1 case and the s-wave shifts towards 
larger doping having its amplitude minimum at the max- 
imum of d-wave. For ^/a = 4 the extended s-wave domi- 



nates and the d-wave either becomes a minor component 
or does not appear at all. Furthermore, in the optimal 
doping whichever symmetry dominates causes the ampli- 
tude of the other minimum i. e, the dominant symmetry 
always expels the other one at the optimum doping. 
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FIG. 4. Same as that in figure 3 except for 9 — i.e 
dx2_y2 -I- Sxy phase that preserves the time reversal symme- 
try. The notable difference is that the usual dx2_y2 and Sxy 
components can mix with each other freely in absence of or- 
thorhombocity, signifying the importance of the phase of 
the non-d-wave symmetry, in contrast to figure 3(a). 



Following the above discussion it is obvious that the 
Fig. 4 represents an exceptional case. Fig. 4 represents 
phase diagram of superconductors having mixed phase 
symmetry like A^^, .{O)^^ + e^'A,^jO)/^-(efc) 
with ^ = (the case oi = 7r/2 is discussed in Fig. 3 
and should be contrasted with Fig. 4). The phase dia- 
gram comprises the amplitudes of the respective sym- 
metry channels as a function band filling p. In striking 
contrast to all the figures Fig.l, Fig. 2 and Fig. 3, there is 
strong mixing of dx2_y2 with Sxy for ^/a =1,3 &4. 
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FIG. 5. Temperature dependencies of the superconducting 
gap in the dx2_y2 and ds^^ channel for their real and com- 
plex mixing for different band fillings (a) p — 0.75 and (b) 
p = 0.9. When the Sxy component has larger Tc, its thermal 
growth is suppressed at the onset of the dx2_y2 component 
(cf. (a)) but that of the dj,2_y2 amplitude is not influenced 
by the corresponding onset of the s^y (cf (b)). In general, 
for 6^ = the gaps open up at a faster rate with decreasing 
temperature than that for 6 — ■k/2. 



In fact mixing between the two symmetries is so strong 
that it is difficuh to find out the predominant symme- 
try for the cases ^/a = 1&3. In this mixed symmetry, 
for e = 7r/2 and ^/a = 4 (cf. Fig. 3(d)), the d-wave 
amplitude is practically zero whereas for 6 — (cf. Fig. 
4(d)) it has strong mixing regime. This is the only mixed 
phase where both of the symmetries at optimal doping 
has large values (see Figs 4(a), (c)) unlike those in figures 
1 to figures 3. The results of this figure thus convincingly 
points out the role of the phase between the two mixing 
symmetries. All the experimentally observed properties 
of cuprates will be consistent with the scenario of Fig. 4, 
including the sign change of the order parameter as well 
as gap nodes. The strong interplay between the two or- 
der parameters of mixed d — s^y symmetry has also been 
reflected in their thermal behaviors (cf. Fig. 5). In Fig- 
ures 5 and 6 we display the temperature dependencies 
of the amplitudes (in eV) of different symmetry order 



parameters for ^/a = 3 as maximum mixing is found in 
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FIG. 6. Temperature dependencies of the superconducting 
gap in the dx2_y2 and 8^2 _y2 channel for their real and com- 
plex mixing for different band fillings (a) p — 0.75 and (b) 
p — 0.9. Unlike the Sxy case (cf. figure 5.), the temperature 
dependencies of the gap amplitudes in the respective chan- 
nels do not influence each other. Similar to that in figure 5, 
for = the gaps open up at a faster rate with decreasing 
temperature than that for 6 = ■k/2. 

this case. When the Sxy component determines the bulk 
Tc, {e.g, at p = 0.75 in Fig. 5(a)) the amplitude of the 
Sxy component is suppressed with the onset of the d-wave 
component. However, when the bulk Tc is determined by 
the d-wave, the amplitude of the d-wave is not affected by 
the onset of the s^y component. This behavior is indeed 
new. In a study of mixed phase with usual d + is phase 



with s as isotropic s-wave, it was shown earlier m9 21 
that the d-wave component gets suppressed with the on- 
set of s-wave but not the reverse. In contrast to Fig. 5, 
the temperature dependencies of the amplitudes of the d 
and 



^x^-{-y 



2 symmetries remain unaffected by each other 
as displayed in figure 6. In general, however, the growth 
of the amplitudes of different symmetries with lowering 
in temperature is faster in case oi 9 — than that for 
6 = 7r/2. This once again emphasize the role of the phase 
0. Temperature dependencies for other values of ^/a is 



qualitatively same as those shown in figures 5 and 6. 
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FIG. 7. Momentum anisotropy of the higher anisotropic 
d-wave symmetry. This higher anisotropic dx2_y2 symme- 
try originates from the fourth neighbour attraction in an 
anisotropic lattice (cf. V4 terms in Eq. (til))- Remarkable 
difference in the fc-anisotropy of this d-wave symmetry com- 
pared to the usual d-wave symmetry is worth noticing. This 
d-wave has 2A{k),nax/kBTc = 5 at p = 0.8. 



So far we have discussed the interplay of order param- 
eters in mixed phases like A^ ^ ^ -I- e* Sa, a = x^ + j/'^ 
or xy. This excluded discussion of some other exotic d 
and Sj;2^y2 symmetries that can arise from the 4*'' neigh- 
bour attraction as discussed earlier in the context of Eqs. 
(|l|,|l). More specifically, a combination of (8*'* + 9*'') and 
(10*'* -|- 11*'*) terms of Eq. (|l|) can give rise to mixed 
pairing symmetries such as, A(fc) = Ad 2_ 2(0)i^''(fc) -I- 
e**A^^,^^,(0)G'"(fc) where F'^{k) = f{k)[l + f'^-y{k) + 
r-«(fc)],V(fc) = g{k)[f-y{k)+J'-y{k) - 1]. These ex- 
otic symmetries are not discussed in the literature. Fol- 
lowing the same procedure as deriving Eq. rtq) one can 
find the gap equation for the components Ad ^ _ ^ (0) and 
As 2 2 (0)j although bit complicated arrives at the same 
gap equation as Eq. (ra) with the pair vertex Vj -^ Vj/2 
and fl — F''-{k), /| — G^{k). Solving the gap equations 
together with the number equation (o) simultaneously 
no mixing between these unconventional d and s wave 
symmetries was found. Within the same parameter as in 
earlier figures (ie, V1/V2 = 0.71), d-wave remains very 
strong at lower dopings (within the range 1 > p > 0.70) 
whereas the s-wave amplitude appears very close to zero 
band filling. Therefore, in Fig. 7 we present the momen- 
tum anisotropy of the unconventional d-wave gap origi- 



nated from 4*'' neighbour attraction. It is clear that gap 
anisotropy is undoubtedly very different form the usual 
nearest-neighbour d-wave symmetry, although basic fea- 
tures of change in sign, nodes etc. remains same as that of 
the ordinary d-wave. This gap symmetry at p = 0.8 gives 
rise to a BCS gap ratio 2A{k)max/kBTc = 5.0 against 
4.29 in case of usual d-wave. Such higher anisotropic d 
wave symmetries will have advantage of avoiding elec- 
tronic repulsion in strongly correlated system like the 
cuprates. 



IV. CONCLUSIONS 

We have studied the superconducting phase with two 
component order parameter scenario, such as, dx2-v2 + 
e'^Sc, where a = xy, x^ +2/^- We showed, that in ab- 
sence of orthorhombocity, the usual dx2_y2 does not mix 
with usual Sx2+y2 symmetry gap in an anisotropic band 
structure. But the Sxy symmetry does mix with the usual 
d-wave for Q — Q. Even in absence of orthorhomboc- 
ity, the higher anisotropic d-wave symmetry mixes with 
higher anisotropic extended s wave symmetry. This is 
obtained by considering longer ranged two-body attrac- 
tive potential in the spirit of tight binding lattice than 
the usual nearest neighbour. This study revealed that 
the dominant pairing symmetry changes drastically from 
d to s like as the attractive pair potential is obtained from 
longer ranged attraction - if the interaction is sufficiently 
short ranged that can be mapped into a nearest neigh- 
bour potential, at low doping, the system is described 
by pure d^2_y2 order parameter. Such consideration of 
longer range attraction has also been revealed by recent 
ARPES data |Q . The role of longer range pair poten- 
tial on pairing symmetry within weak coupling theory of 
superconductivity has thus been established. We showed 
that the momentum distribution of the higher anisotropic 
d-wave symmetries is quite different from the usual d- 
wave symmetries. We found that the typical interplay 
in the temperature dependencies of these higher order 
d and s wave pairing symmetries can be different from 
what is known. In brief, we believe such study of higher 
anistropic symmetries is potentially important and will 
stimulate further studies in contrast to the usual d and s 
wave symmetries. 
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